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Abstract 

We consider the following dynamic Boolean model introduced by van den Berg, Meester and 
White (1997). At time 0, let the nodes of the graph be a Poisson point process in M. d with 
constant intensity and let each node move independently according to Brownian motion. At 
any time t, we put an edge between every pair of nodes whose distance is at most r. We study 
three fundamental problems in this model: detection (the time until a target point — fixed or 
moving — is within distance r of some node of the graph); coverage (the time until all points 
inside a finite box are detected by the graph); and percolation (the time until a given node 
belongs to the infinite connected component of the graph). We obtain precise asymptotics for 
these quantities by combining ideas from stochastic geometry, coupling and multi-scale analysis. 

Keywords and phrases. Poisson point process, Brownian motion, coupling, Minkowski dimen- 
sion. 
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1 Introduction 

In a random geometric graph, nodes are distributed according to a Poisson point process in R rf of 
intensity A (i.e., nodes are uniformly distributed, with A nodes in expectation per unit volume), 
and an edge is placed between all pairs of nodes whose distance is at most r. These graphs, also 
known as "continuum percolation" or the "Boolean model," have been used extensively as models 
for communication networks [51 114} [T5]. While simple, these models capture important qualitative 
features of real networks, such as the phase transition in connectivity as the intensity A of the 
Poisson point process is increased. 

In many applications, the nodes of the network are not fixed in space but mobile. It is natural 
to model the movement of the nodes by independent Brownian motions [3JEJHI]; we call this the 
mobile geometric graph model. We focus on three fundamental problems in this model: detection 
(the time until a target point — which may be fixed or moving — comes within distance r of a node 
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of the graph); coverage (the time until all points in a large set are detected); and percolation (the 
time until a given node is connected to the infinite connected component). 

We now give a precise definition of the model and state our results. Motivation and related work 
are addressed in Section [2j 

The Mobile Geometric Graph model 

Let Ho = {Xi}i be a Poisson point process on R rf of intensity A. To avoid ambiguity, we refer to the 
points of a point process as nodes. We let each node Xj move according to a standard Brownian 
motion {d(s)) s >o independently of the other nodes, and set H s = {Xi + Ci( s )}i to be the point 
process obtained after the nodes of Ho have moved for time s. By standard arguments [3] it follows 
that H s is again a Poisson point process of the same intensity A. 

At any given time s we construct a graph G s by putting an edge between any two nodes of Tl s 
that are at distance at most r. In what follows we take r to be an arbitrary but fixed constant. 
There exists a critical intensity A c = \ c (d) such that if A > A c , then a.s. there exists a unique 
infinite connected component in G s , which we denote by Coo(s), while if A < A c then all connected 
components are finite a.s. |14| 119] . 

Statement of results 

Detection. Consider a "target" particle u which is initially placed at the origin, and whose position 
at time s is given by g(s), which is a continuous process in R d . We are interested in the time it 
takes for u to be detected by the mobile geometric graph, namely how long it takes until a node of 
Ilo has come within distance at most r from u. More formally we define 

T dct = inf {* > : g(t) G B(Xi + Ci(*),r)}, 

i 

where the union is taken over all the nodes {Xi} of Ilo and B(x,r) denotes the ball of radius r 
centered at x. In Section [3] we prove the following theorem, which extends previous classical results 
on detection for non- mobile particles u \24:\ [9] . 

Theorem 1.1. In two dimensions, for any fixed A and any g independent of the motions of the 
nodes of Ho, we have that 

P [T dct >t}< exp f-27rA^(l + o(l))) . 

In addition, if g is an independent Brownian motion then the above bound is tight, i.e., 

P [T det > t] = exp ^-2ttA^(1 + o(l))\ . (1) 

Remark 1.2. Theorem II. II is stated for d = 2. In Section [3] it is extended to all dimensions d > 1 
(see Theorem 13. 5p . where the tail of T det is shown to be exp(— @{yt)) for d = 1 and exp(— 0(i)) 
for d > 3. All these results exploit a connection between the detection time and the volume of the 
Wiener sausage, as explained in Section [3j 

Remark 1.3. From |24^ [9], we have that for d = 2 the result in ([1]) also holds when u does not 
move (i.e., g = 0). Thus Theorem 11.11 establishes that, asymptotically, the best strategy for a 
particle u (that is not informed of the motion of the nodes of Ilo) to avoid detection is to stay 
put. In Section [3] we show that this is true for d = 1 for any fixed t (not only asymptotically) and 
conjecture that this holds in higher dimensions as well (see Conjecture I3.6p . 
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Coverage. Let A be a subset of R d . We are interested in the time it takes for all the points of A 
to be detected. We thus define 

T cov (,4) = inf {t > : A C \J\J B(Xi + Q(s),r)}. 

s<t i 

For R E R_|_, let be the cube in M. d of side length R. A natural question proposed by Konstan- 
topoulos [H] is to determine the asymptotics of E [T cov (Qjt)] as R — > oo. 

In Section [5] we prove Theorem 11.41 below, which gives the asymptotics for the expected time to 
cover the set Qr as R — > oo and shows that T cov (Q#) is concentrated around its expectation. We 
write / ~ g as x — > oo, when — > 1 as x — > oo. 

Theorem 1.4. We have that as R — > oo 



ET C0V (Q R ) ~ { 



^(log^ /o rd=l 

^ log i? log log ford = 2 and r 1 in Probability, 

/ord > 3 E[T cov (g K )] 



where c(d) = d and V stands for the Gamma function. 

Remark 1.5. Instead of covering a whole cube, we could ask for the coverage time of other sets. 
We prove Theorem 11.41 in this general setting in Section HJ for instance, we show that E [T cov ] 
for a line segment of length R is smaller than E[T cov ((5r)] by a factor of 1+ °} l " > and also obtain 
asymptotics for fractal sets (see Theorem I4.2p . 



Percolation. Let u be an extra node initially at the origin and which moves independently of the 
nodes of IIo according to some function g. We now investigate the time it takes until u belongs 
to the infinite connected component. We denote this time by T perc , which can be more formally 
written as 

T pcrc = inf{t > : 3y G C^t) s.t. \\g{t) - y\\ 2 < r}. 

The detection time clearly provides a lower bound on the percolation time, so we may deduce from 
Theorem 1 1.1 1 and Remark 1 1 . 2 1 ab ove that P [T perc > t] is at least exp (— 0(t/ log t)) for d = 2 and at 
least exp (— 0(t)) for d > 3, when u is non-mobile or moves according to an independent Brownian 
motion. We will prove the following stretched exponential upper bound in all dimensions d > 2 in 
Section [5j 

Theorem 1.6. For all dimensions d>2, if X > X c (d) then there exist constants c and to, depending 
only on d, such that 

P [Tpcrc > t] < exp (y~ c lo J+6/d^ j ' f° r al1 1 - f °- 

This holds when u is non-mobile or moves according to an independent Brownian motion. 

Theorem 11.61 is the main technical contribution of the paper; we briefly mention some of the ideas 
used in the proof. The key technical challenge is the dependency of the G s 's over time. To overcome 
this, we partition Mr into subregions of suitable size and show via a multi-scale argument that all 
such subregions contain sufficiently many nodes for a large fraction of the time steps. This is the 
content of Proposition 15.21 which we believe is of independent interest. This result allows us to 
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couple the evolution of the nodes in each subregion with those of a fresh Poisson point process of 
slightly smaller intensity A' < A which is still larger than the critical value A c . After a number 
of steps A that depends on the size of the subregion, we are able to guarantee that the coupled 
processes match up almost completely. As a result, we can conclude that there are Q(t/A) time 
steps for which the mobile geometric graph contains an independent Poisson point process with 
intensity A" > A c . This fact, which we believe is of wider applicability, is formally stated in 
Proposition 15.11 This independence is sufficient to complete the proof. 

Finally, to illustrate a sample application of Theorem 11.61 we consider the time taken to broadcast 
a message in a network of finite size. Consider a mobile geometric graph in a cube of volume n/A 
(so the expected^ number of nodes is n). Since the volume is finite, we need to modify the motion 
of the nodes to take account of boundary effects: following standard practice, we do this by turning 
the cube into a torus (so that nodes "wrap around" when they reach the boundaries). Suppose 
a message originates at an arbitrary node at time 0, and at each integer time step t each node 
that has already received the message broadcasts it to all nodes in the same connected component. 
(Here we are making the reasonable assumption that the speed of transmission is much faster than 
the motion of the nodes, so that messages can travel throughout a connected component before it 
is altered by the motion.) Let T^ c denote the time until all nodes have received the message. We 
prove the following result in Section [6j 

Corollary 1.7. In a mobile geometric graph on the torus of volume n/A with any fixed A > A c , the 
broadcast time T^ c is 0(log n(log log n) 3+6 / rf ) w.h.p. in any dimension d>2. 

2 Motivation and related work 

A possible application of our results is in the study of mobile ad hoc networks, where nodes moving 
in space cooperate to relay packets on behalf of other nodes without any centralized infrastructure. 
This is the case, for example, in vehicular networks (where sensors are attached to cars, buses or 
taxis), surveillance and disaster recovery applications (where mobile sensors are used to survey an 
area), and pocket-switched networks based on mobile communication devices such as cellphones. 

Random geometric graphs have long been used as a model for static wireless networks, and by now 
their structural and algorithmic properties are rather well understood mathematically. We refer 
the reader to the book [19] for extensive background on random geometric graphs and to I22| for a 
selective survey of applications. We should mention that in many papers the model is defined over 
a torus 5 C l rf of finite volume n/A, so that the expected number of nodes is n. We choose to 
work in the infinite volume which is mathematically cleaner, but most results obtained there 
can be adapted to finite volumes with a little technical work; see Section [6] for an example. In 
the finite setting it makes sense to talk about the random geometric graph being connected, which 
is an important property when the nodes are static. This occurs at a sharp threshold value of A, 
which however grows with n [TJ \T7\ [18] . 

The scope of mathematically rigorous work with mobile nodes is much more limited, and there is 
as yet no widespread agreement on an appropriate model for node mobility. The model we use 
in this paper is equivalent to the "dynamic boolean model" introduced by Van den Berg et al. [3] 
(who also proved that almost surely an infinite component exists at all times if A > A c ). We point 
out that in this model, in contrast to many others, node mobility is fixed and does not depend on n 
(the number of nodes in a finite network). 

The result can be adapted to the case of a fixed number of nodes n using standard "de-Poissonization" argu- 
ments |19| . See the Remark following the proof of Corollary 11.71 in Section [6] 
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The detection time was addressed by Liu et al. [12], assuming that each node moves continuously 
in a fixed randomly chosen direction; they show that the time it takes for the network to detect 
a target is exponentially distributed with expectation depending on the intensity A. Also, for 
the special case of a stationary target, as observed in [9j [11] the detection time can be deduced 
from classical results on continuum percolation: namely, in this case it follows from [23] that 
Pr [Td c t >t]= exp(— A E [vol (W r (t))]), where W r (t) is the "Wiener sausage" of radius r up to time t. 
This volume in turn is known quite precisely [23^12]. In a different model, in which nodes perform 
continuous-time random walks on the square lattice, Moreau et al. |15j (see also [6j) establish that 
the best strategy for a target to avoid detection is to stay put. (In these papers, a target is said to 
be detected when there exists a node of the graph at the same lattice location as the target.) The 
analysis of the coverage time was suggested as an open problem by Konstantopoulos |llj . 

The percolation time was first studied in [22j . which also examined detection for moving targets. 
This paper is a strengthened version of [22], with tighter results on percolation and detection as 
well as the addition of coverage. 

The question of broadcasting in mobile graphs has been studied by several authors in a setting 
where messages travel only to immediate neighbors at each time step. Clementi et al. [5] establish 
tight bounds for the broadcast time in this setting assuming that either the intensity A or the 
range of motion of the nodes grows with n. The case of smaller intensities and bounded range of 
motion was studied by Pettarin et al. [21]. A problem similar to broadcast was studied by Kesten 
and Sidoravicius [ID], who derived the rate at which an infection spreads through nodes that are 
performing continuous-time random walks on the square lattice. 

3 Detection time 

In this section we give the proof of Theorem 11.11 We first state a generalization of a well-known 
result [23], which we will use in several proofs; we include its proof here for the sake of completeness. 

Lemma 3.1. Suppose that u starts from the origin at time and its position at time s is given by 
a deterministic function g(s). Let W g (t) = L) s <tB(g(s) — £(s),r) be the so-called "Wiener sausage 
with drift" up to time t. Then, for any dimension d > 1, the detection probability satisfies 



Since the (j's are independent we have that $ is a thinned Poisson point process with intensity 
given by 




Proof. Let $ be the set of points of IIo that have detected u by time t, that is 

$ = {Xi G n : 3s < t s.t. g(s) G B{X { + &(s),r)}. 



A(x) = A P [x G U s < t B(g(s) - C(s), r)] , 
where £ is a standard Brownian motion. 

So for the probability that the detection time is greater than t we have that 




P [Tdet >t} = exp(-A / P [x G U s < t B(g(s) - ((s),r)} dx) 



exp(-AE [vol(U s <t J B(< 7 ( S ) - C00,O)D = exp(-AE [vol {W g {t))}). 



□ 
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Remark 3.2. The preceding lemma implies that when the motion g of u is random and independent 
of the motions of the nodes of the Poisson point process IIo then 

P [Tdet > t] = E [exp(-AE [vol(W s (t)) | (g(s)) s < t })} . 

Remark 3.3. We note that the above proof can be easily generalized to show that the time T^ et (K) 
until we detect some point in a compact set K C R rf satisfies 

P [T dct (K) >t} = exp(-AE [vol (U s < t (K r - C(«)))D, (2) 

where K r stands for the r-enlargement of K, i.e. K r = L) x ^KB(x,r). 

From Lemma[3J]we see that estimating P [T dct > t] translates to deriving estimates for E [vol (W g (t))]. 
When u does not move (i.e., g = 0), it is well known that in two dimensions [23} 12] 



E [vol (W (t))] 



2vrf 
log* 



(l + o(l)). 



The following Lemma implies that E [vol (W s (i))] > E [vol (Wo(i))] (1 — o(l)) for any deterministic 
continuous g. 

Lemma 3.4. Let Q be a standard Brownian motion in two dimensions and let g be a deterministic 
continuous function, g : M + — > IR 2 . Let W g (t) = U s <tB(g(s) — ((s),r) be a Wiener sausage with 
drift g up to time t. We then have that as t — > oo 

E[vol(W g (t))]> gl(l- (l)). 

Proof. We may write 

E [vol(W a (t))] = / P [y G U s < t B(g(s) - ((s),r)] dy = [ P [r B(tf , r) < t] dy, 
where ta is the first hitting time of the set A by g — (. Define 

Z y = f l(g(s)-((s)e B(y,r))ds, 
Jo 

i.e., the time that the process g — C spends in the ball B(y,r) before time t. It is clear by the 
continuity of g — ( that {Z y > 0} = {t b ^ ^ < t}. Clearly P [Z y > 0] = e[z | z >o] 



and for the first 



moment we have 



E[Z y }= [ P\g(s)-((s)eB(y,r)] ds 
Jo 



1 IN-aWllj 
e 2 S dzds 

2<KS 



JB(0,r) 



1 |z+y-g(s)ll2 

e 2 S dzds. 

2tts 



10 JB(y,r) 

For the conditional expectation E [Z y \ Z y > 0] , if we write r for the first time before time t that 
g — C hits the boundary of the ball B(y,r), denoted by dB(y,r), then we get 

E[Z y \Z y >0) 



E 

(t-r)Vl 



pt—T 

/ 1(50 + t) -C(s + t) G B(y,r))ds 
Jo 



< 1 + E 

< 1 + max 

< 1 + 



+ t) - ((s + r) G r)) 



E 



xedB{y,r) J 1 J B(y,r) 
t 



-g(s + T)-3H-g(T)||g 



2vrs 



(is 



1 yS(s/,r) 



J-dzds<l + r 2 ^. 
2vrs ~ 2 
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So, putting everything together we obtain that 

E [ZJ JB (° 



E[vo 



Kw 9 (t))} = [ 



and hence as t — > oo 



E [Zy \Zy>0] 



dy > 



I II zJ rV~ Q( s ) II 2 

Jo Ib(o,t) Is? 2Ts e Ta d y dzds 



2ntr 2 



1 + r 



2 log* 



2 + r 2 log t 



2vrt 



E[vol(W g (t))]> —(l- (l)). 



□ 



We are now ready to prove Theorem I l.li 

Proof of Theorem From Remark 13.21 we have 

P [Tdet > t] = E [exp(-AE[vol(W fl (f)) | (g(a))«<t])] 
where 5 is independent of £. By Lemma |3,4| we have the upper bound 



P [T d et > t] < exp ^-2vrA— (1 - o(l)) 1 , as t -> 00. 

So it remains to show the lower bound on this probability for the case when g is a standard Brownian 
motion independent of the motions of the nodes of Ilo- Letting R = logi, it is clear that 

P [Tdet > t] > P [u stays in B{0, R) for all s < t, T B(0>R) > t] , (3) 

where Tb(q,r) is the detection time of the ball B(0, R), i.e., 

T B (o,R) = mf{s > : 3i s.t. B(X { + &(s),r) D B(0,i2) / 0}. 

Since the motions of u and the nodes of Llo are independent, we get that 

P [u stays in B(0, R) for all s<t, T B{0 R) > t] = P [u stays in B(0, R) for all s < t] P [T B{0 R) > t] 

(4) 

From Remark 13.31 with K = B(0, R) we get 

P [T B (o,R) >t)= exp(- A E [vol (U s < t B{0, R + r)- £(*))]) 

and writing U s < t (I?(0, R + r) — ((s)) = R U s < t B (—^7^, 1 + jjj and for i? large enough we get 
that, for all s, 
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Of) 

R 



1 Cfl 



,2 • 



For any i>Owe have by Brownian scaling that 

COO 



E 



vol MJ B 



,s<t 




E 



vol 



U b(c(s% 



8'<4t 

— a 3 



where £ is a standard Brownian motion. So finally, using the asymptotic expression for the expected 
volume of the Wiener sausage in two dimensions [21 [23], i.e., E [vol (U s <t-B(C(s), x))] ~ ^u^l as 
t — > 00, for any x independent of t, we get that 



E [vol (U 8 < t B(0, R + r)- £(«))] ~ R 2 2n 



t/R 2 

loET f — log R 2 



2vr- 



logt' 



Hence, 



P [T B (o,R) >t)= exp ^-2vrA^(l + o(l))J , as t -> oo. 



Thus we only need to lower bound the probability that u stays in the ball B(0,R) for all times 
s<t. 

For any t > R 2 and any dimension d > 1, we have by [I] that 

P [u stays in B(0, R) for all s < t] > exp(-ct/i? 2 ), (5) 
for a positive constant c and hence, since i? = log t, we get that 

P [T det > t] > exp ( "2vrA i -^(l + (1))) exp (^7^2) = exp (-2ttA^(1 + o(l)) 

□ 



General dimensions and the Wiener sausage 

For d = 1, the volume of W g (t) can be computed from the maximum and minimum values of g — £ 
via the formula 

vol (W ff (i)) = 2r + max(<?(s) - C(s)) - miri( ff (s) - C(a)). 

Let i* and be the random times in the interval [0, t] at which — £ achieves its maximum and 
minimum values respectively. Then we have 

E [vol (W g (t))] > 2r + E [(<#*) - C(t*)) - (<?(**) - C(**))] 

= 2r + E [-C(t*) + C(Q] = E [vol (W (t))] , (6) 

where E [<?(£*)] = E [<?(£*)] holds since t* and i* have the same distribution. Thus, for d = 1 the 
inequality E [vol (W g (t))] > E [vol (Wo(£))] holds for all fixed t; for d = 2 Lemma I3T41 gives this 
inequality only asymptotically as t — > oo. 

For dimensions d > 3, the proof of Lemma 13.41 can be used to obtain the following weaker result: 
there exists a positive constant c such that 

E [vol (W g {t))} > cE [vol (W (t))] . (7) 
The expected volume of the Wiener sausage with g = is known to satisfy |23t [2] 

^ + 2r for d = 1 



V (t) =E [vol (W (t))] 



g(l + o(l)) ford = 2 (8) 



I, c(d)r d " 2 i(l + o(l)) for d > 3, 

where c(d) = ^572 • (For d = 1 the quantity above follows from well-known results for Brow- 
nian motion [16]; namely E [max s < ( C( s )] = — E [min s < t £(s)] = y~-) Hence plugging (|6HH|) into 



Lemma 13.11 gives an upper bound for P [Tdet > *] when d = 1 and d > 3. Regarding the lower 
bound, when g is an independent Brownian motion, the same strategy as in the proof of Theo- 
rem [1J] works for d = 1 and d > 3 provided we set R properly. For d = 1 it suffices to take R = t 1//3 ; 
for d > 3 we can set R = r. Then we obtain a positive constant c\ such that, as t — > 00, 



P[T B{0R) >t]>\ exp(-^f(l +o( l))j ford=l 
[exp (- Cl r d - 2 t(l + o(l))) for d > 3. 

This together with (JH) and ([5]) give us the following theorem, which holds in all dimensions d > 1. 
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Theorem 3.5. Let g be a continuous process in M. d , for d > 1. If g is a deterministic continuous 
function, then we have 

-i log P [Tdet > t] = E [vol (W^t))] . 

If g is random but independent of the motion of the nodes o/Hq, we obtain a positive constant (5\ 
such that as t — > oo 



1 



V (t) ford=\ 

\ log P [T^ >t}>< (1 - o(l))V (t) ford = 2 

[1 - o(l))p x V Q {t) ford>3, 

where Vo(t) is defined in (jSJ). 

If g is a standard Brownian motion, then we obtain a positive constant fa such that as t — >■ oo 

"A l0gP[Tdct>t] -\ (l + o(l))AV (t) ford>3. 

We conclude this section with the following conjecture. 

Conjecture 3.6. For d > 2, any /fared t, and any function g, we have that 

E[vol(W 5 (t))] >E [vol (Wo (*))]■ 

Note that Theorem 13.51 establishes the conjecture for d = 1, and establishes it asymptotically as 
t — > oo for d = 2 provided g is continuous. A related statement for random walks was proved in 
[15] : see also [U Corollary 2.1]. 



4 Coverage time 

In this section we will prove a more general version of Theorem Ol Let A be a subset of R d . For 
i? £ M + we define the set iL4 = {i?a : a £ A}. We recall the definition of Minkowski dimension, 
which can be found, e.g., in [13]. 



Definition 4.1. Let A be a non-empty bounded subset of Mr. For e > let M(A,e) be the smallest 
number of balls of radius e needed to cover A: 

M(A, e) = min jfc : A C (J B(xi, e) for some x-i G M d | . 

The Minkowski dimension of A is defined as 

log M(A,e) 



dim a/ (^4) = lim 



e^o log e" 
whenever this limit exists. 



We now proceed to state the more general version of Theorem 11.41 



9 



Theorem 4.2. Let A be a bounded subset of M. d of Minkowski dimension a. We have that as 
R — > oo 



ET cov (RA) ~ { 



f^(logi?) 2 ford=\ 



T ( RA) 

2^ log R log log R for d = 2 and ™ v { ' ->• 1 in probability, 

/ord > 3 E[r cov (iL4)] 



where c(d) = d and V stands for the Gamma function. 

2-kT. 

Proof. In the proof we will drop the dependence on RA from T cov (RA) and E [T cov (iL4)] to simplify 
the notation. We will carry the proof for the case d = 2 only and discuss how to adapt the proof 
for other dimensions at the end. 

Let M(A, e) = min{fc > 1 : 3 B±, . . . ,B^ balls of radius e covering A}; then it is easy to see that 
M(RA, e) = M(A, -^). By the assumption that A has Minkowski dimension a, for any S > we 
can find eo small enough such that 

e ~a+8 < M ^ £ ) < -a-S^ for any £ < ^ (g) 

We will first show that 

hmsup-j— 51 — 1 B - L 

/?^oo ^ log # log log # 

To do so, we are going to cover the set RA by M = M(RA, e) balls of radius < e < r. 

From (J9]) we get that 

/ g \ — a+5 / £ \ —a— 5 

(-^J < M < i^—j for R sufficiently large. (10) 

Let Zt be the number of balls not covered by the nodes at time t. It is clear that {T cov > t} C 
{Zt > 1}. For the first moment of we have 

E [Zj] < Af P [a given ball B(x, e) is not covered by time t] . 

The probability that a ball B(x,e) is covered by time i is lower bounded by the probability that 
a node of the Poisson point process Ho has entered the ball B(x,r — e) before time t. Hence, 
P [B(x, e) is not covered by time t] is at most the probability that x has not been detected by time 
t by a mobile geometric graph with radius r — e. From Lemma 13.11 we obtain 



P [B(x, e) is not covered by time t] < e 



-AE[vol(WV_ e (t))] 



where W z , p (t) = U s < t B(z + ((s), p). 

We are now prove the upper bound. Let 5' > be small. For t large enough we have that (see 
©) 

(1 _ S ')2w J- < E [vol(Wb,r- e (<))] < (1 + S')2ir-L (11) 
log t log t 

and hence, 

B[Z t ] < A/e~ 2wA(1 ~ 5,) i5k (12) 
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By Markov's inequality we have that P [T cov > t] < E [ZJ. Also, 

E[T cov ]= / P[T cov >i] dt<t*(R)+ / Me^^ 11 ^'^ dt, (13) 



where t*(R) satisfies 



Mexp f -2vrA(l - <^) - **^L ) = 1 and U(R) > e, for sufficiently large R. (14) 
\ log U (R) J 

We claim that the last integral appearing in (fl~3j) is o(t*(R)). To see this set c = 2-7rA(l — 5') and 
use a change of variable, x = ^|^, which gives £j| > 2 io g t — 4 io g S ^ or * l ar § e enough. Hence setting 



t (R) 

log t [fe ^ ne integral is upper bounded by 



POO 

/ 4A/(logx)e" c:E (ix < c'(l/x* + logx*) = o{t*(R)). 
So we finally obtain that 

E [T cov ] < t»(i2)(l + o(l)), as R -> 00. 

From dH]) and (JTUJ) we get 

-R i (-^) 

(a — (5) log — < c- -— < (a + 5) log — for i? large enough, 

e log t* (R) e 

and thus we conclude that 

hmsup - 1n _ BW1n _ B < 1, 

which follows by letting 5 and 5' go to 0. 
We now proceed to show the lower bound 

liminf „ , E £H, > L 

sfxiog-ffiogiog-R 

To do so, we are going to use the equivalent definition of Minkowski dimension involving packings 
[TBI Chapter 5]. Letting 

K(A, e) = maxjfc > 1 : 3 B\, ■ ■ ■ , disjoint balls of radius e centered in ^4}, 

it is clear that K(RA, e) = K(A, -^). For 5 > there exist K = K(RA, 1) disjoint balls with centers 
in RA and radius 1 satisfying 

R a-8 < K < R a+S^ for R large enough 

So, for R large enough, we can pack the set RA with points xi, ■ ■ ■ , xk (the centers of the balls) 
that are at distance at least 2 from each other. Let Ut denote the number of centers x±, ■ ■ ■ ,xk 
that have not been detected by time t. Obviously we have that {T cov > t} D {Ut > 1}. 

Recall that the Wiener sausage W z>r (t) = L) s <tB(z + (,(s),r) in two dimensions satisfies, for 5' > 0, 

(1 - 5')2ir-^— < E [vol (W z r (t))\ < (1 + 5')2tt—, for t large enough. (15) 
log t log t 
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Let e > be small and let t* = t*(R) > e satisfy the equation 



t* OL € S , _ . . 

77 log R. (16) 



logt* 2vrA(l + (5 / ) 
Applying the second moment method to the random variable Ut* we obtain 



P [T cov > f] > 



so in order to obtain a lower bound for P [T cov > t*] it suffices to lower bound the first moment of 
Ut* and upper bound its second moment. We will show that P [T cov > t*] > jqr^ni ; hence we will 
get that E[T cov ] > t*j^. 

We have that E [Ut*] = Yli=i P I x « n °t detected by time £*], and using Lemma [37X1 we obtain that 



P [xi not detected by time t*] = exp(-AE [vol(W Xi )]), 

where W x = W x , r (t*) = U s < t * B{x + ((s), r). 

Obviously E [vol(Wz)] is independent of x, and hence we get that 

E[U t *] = Kexp(-AE[vol(VF )]). (17) 

Now, for the second moment of Ut* we have 

K 

E [U£] = ^2 ^2 P X J not detected b y time **] + E Pf] ( 18 ) 

i=l j'^i 

and using Remark 13.31 we get that 

P [xi,xj not detected by time t*] = exp(— AE [vo^W^ U W^.)]). 

(Note that the two Wiener sausages W x . and W Xj use the same driving Brownian motion.) 
Writing 

vol(W Xi U W Xj ) = vo\(W Xi ) + vol(W Xj ) - vol(W Xi n W Xj ), 
equation (fTKj) becomes 

N 

E [E/£] < exp (-2AE [vol(Wb)]) ^ ^exp(AE [vol(W Xi n W Xj )] ) + E [U t *] . (19) 

i=l j^i 

Thus it remains to upper bound E [vol(W Xi D Wj.)] for all i and j. If \\x{ — xj\\2 < (logi?) 2 , then 
we may use the bound vol(W Xi H W^) < vol(Wa-). 

Recall from (|16p that t*(R) = ©(log i? log log it!). The idea is that if X\ and xj are at distance 
greater than (logi?) 2 apart, then it is very unlikely that the two sets W Xi and W x - will intersect. 
Specifically, when ||sGj — Xj\\2 > (logR) 2 it is easy to see that the probability that the two sausages, 
W Xi and W x -, intersect is smaller than the probability that a 2-dimensional Brownian motion has 
traveled distance greater than ^(logi?) 2 in t* time steps, and this last probability is bounded above 
by ce -c( - logR ) 2 by the standard bound for the tail of a Gaussian. 
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When ||sEj — Xj\\2 > (logR) 2 , writing Si = B(xi,R) for the ball of radius R centered at Xi and 
defining inductively S k = B(x i ,2 k - 1 R) \ B( Xi ,2 k - 2 R) for all k > 2, we can split the volume of 
W Xi n W x . as follows: 



E [vo\{W Xi nW Xj ))=f^-E [yol(W Xi D W Xj n S n )] 

n=l 

oo 

< crfe^*) 2 +J2c'2 2n R 2 e- c ' 22nR < c"iT M 



n=2 

where c, c', c" and M are all positive constants. The first part of the first inequality follows from 
the discussion above, namely that if the intersection is nonempty, then the Brownian motion must 
have traveled distance greater than ^(logi?) 2 in less than t* steps. If this has happened, then we 
simply bound the intersection of the two Wiener sausages in the ball B(x,R) by the volume of 
the ball. The second part of the first inequality follows by the same type of argument, since now 
in order to have a nonempty intersection in the set S n , the Brownian motion must have traveled 
distance at least 2 n ~ 2 R in less than t* steps, which again is exponentially small. 

Finally, the sum appearing in (|19p is bounded above by 

fXifllx* " x sh < (iogi?)V E[vol(Wo)1 + EE 1 ^ - x ih > (logi?) V AirM (20) 

i=l jj^i i=l jjti 

By (|15p and the definition of t* given in (|16|) we get that (j20|) is bounded from above by 

Cl K{\ogRfR a - s -' + K 2 e c " XR ~ M 

and hence 

E [Ul] < exp(-2AE [vol(Wo)]) (ci#(log R) 4 R a - s ~ e + K 2 e c " XR ~ M \ + E [U t *\ . (21) 
Therefore, putting all the estimates together we get that 

Pin >nl>^M£> I 

11 J " E[U 2 ] ~ e c"A^ + -^((logi?)4 J R«-5-6 + e AE[vol(W )])- 

Using the lower bound K > R a ~ 5 , the upper bound for the expected volume from (|15f) and the 
definition (fT6l) of t* , we deduce that 



and thus E [T cov ] > pp^rn £*■ Since t* satisfies (fT6|h we deduce that 

liminf E > 1 ~ I ~ I 

fl->oo ^log^Rloglogi? - 1 + 5' 

and hence, letting e, 5 and 5' go to 0, we get that 

liminf , E i?H, > L 
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So, we have shown that 



E[T C( 



a 
2^A 



log R log log R, as i? — > oo. 



(22) 



Now, for d = 2 it only remains to show the last part of the theorem, namely that E j^ ov i converges 
to 1 in probability as R — > oo. For any 7 > we have that 



E[T C( 



> 7 



P [T cov > (1 + 7) E [T cov ]] + P [T cov < (1 - 7) E [T a 

< E [Z (1+7)E[Tcov] ] + P [T cov < (1 - 7) E [T„„ 
From (I12p and the definition of M we have that 

E \Z t ) < R a+5 exp ( -2ttA(1 - <5')t-^ . 

V log*/ 

Plug ging in i — (1 + 7) E [Tcov]) using (1221) and taking 5 1 sufficiently small gives that 



B[Z 



(l+ 7 )E[T cov ]J 



0, as R — > 00. 



For e,S,5' small enough we get that (1 — 7) E [T cov ] <t*, so 

P [T cov < (1 - 7) E [T cov ]] < P [T cov < t*] < P [E/f = 0] < 1 - (E l^*Jf = 



E [U£] 



Hence we get the desired result that 



E[T C , 



>7 



0, as R — > 00. 



For dimensions d 7^ 2, the same arguments carry through by employing the proper expression for 
the expected volume of the Wiener sausage given in ([8]). Then, we need to set t+(R) and t*(R) 
correspondingly. From (fll"]) and (fT6|) . it suffices to set to satisfy 



exp A(l - 5') 



M for d = 1 



7T 



exp ( A(l - 5')c{d)r d - 2 U(R) ) = M for d > 3, 



and t* to satisfy 



A(l + 5 / )\/^^ = (a-e-5)logR for d = 1 



7T 



A(l + (J'Xd)^- 2 *^^) = (a - e - 5) log ii for d > 3. 



□ 

Remark 4.3. While the limit defining Minkowski dimension in Definition 14.11 may not exist, the 
corresponding limsup is denoted by dimjvf(^4) and always exists. The proof of Theorem 14.21 also 
shows that for d = 2 

lrmsup^— — — = dim M (A) 

R ^°° 2^A log log log i? 

and similarly for liminf and other dimensions. 

Remark 4.4. The estimates in the proof of Theorem 14.21 actually imply that a.s. e[t (iM)] ~~ ^ ^ 
as R — > 00. 
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5 Percolation time 



In this section we give the proof of Theorem 11.61 We will observe the process (Gi)«>o in discrete 
time steps i = 0, 1, . . . in order to be able to apply a multi-scale argument. For a nonnegative 
integer i we define the event Ji that u does not belong to the infinite component at time i; more 
formally, 

Ji = {u <£ U y&Coo(i) B(y,r)}. 
Then it is easy to see that, for all t, we have 

P [Tperc > t] < P 



n L * J J,- 

1 1 i=0« y * 



We define Ql to be the cube with side length L centered at the origin and with sides parallel to 
the axes of M. d . We tessellate Ql into subcubes of side length £ < L , which we call cells. We now 
state two key propositions that lie at the heart of our argument. 

The first proposition says that, provided every cell of the tessellation contains sufficiently many 
nodes, then we can couple the positions of these nodes after sufficiently many steps with the nodes 
of an independent Poisson point process of only slightly smaller intensity on a smaller cube. We 
prove this proposition in Section [5.11 

Proposition 5.1. Fix K > I > and consider the cube Qk tessellated into cells of side length 
I. Let #o be an arbitrary point process at time that contains at least f3£ d nodes at each cell of 
the tessellation for some (3 > 0. Let $a be the point process obtained at time A from <I>o after the 
nodes have moved according to standard Brownian motion for time A. Fix e £ (0, 1) and let E be an 
independent Poisson point process with intensity (1 — e)/3. Then there exists a coupling of 'E and $a 
and constants c\,C2,c^ depending only on d such that, if A > an d K' < K — c<i\] A log e _1 > 0, 
then the nodes of 5 are a subset of the nodes of inside the cube Qk 1 with probability at least 

l-^exp(-c 3 6 2 /3f*). 

The second proposition, which we prove in Section 15.21 says that the above condition that each 
cell contains sufficiently many nodes is satisfied at an arbitrary constant fraction of time steps with 
high probability. 

Proposition 5.2. Lett > be a sufficiently large integer and£,e € (0, 1) be two constants. Suppose 
that the cube Ql, for L = t, is tessellated into cells of side length £, where £ d > Clog 3 t for some 
sufficiently large constant C. For i = 0, 1, . . . let 

Ai = {at time i all cells contain > (1 — £)\£ d nodes of Ilj}. 

Then there exists a positive constant c such that 



t-i 



J3l(^)>(l-e)t 



i=0 



Proof of Theorem \1.6[ Let u be a node that is at the origin at time independent of the nodes of 
Ho- We assume that u is non-mobile; the proof can easily be extended to mobile u using translated 
cubes that track the motion of u as in [22^ Section 4]. 
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Let t be an integer sufficiently large. We consider the cube Ql, for L = t. Set Ht to be the event 
that u has never been in the infinite component from time to t — 1. More formally, we define 

H t = nlzlJi = Pilzliu i U yeCoo{l) B(y,r)}. 

We say that a cube Ql has a crossing component at a given time % if among the nodes in Ql there 
exists a connected component that has a path connecting each pair of opposite faces of Ql- (A 
path connects two faces oi Ql if for each face there is at least one node of the path within distance 
r of the face.) We then define Ht to be the event that u has never been within distance r of a 
crossing component of Ql from time to t — 1. Let Kt be the event that, in each step from 
to t — 1, there exists a unique crossing component oi Ql and it intersects the infinite component. 
Therefore, if Kt holds and u belongs to a crossing component of Ql at some time step from to 
t — 1, then at the same time step u will also belong to the infinite component. We can then conclude 
that K t n Ht C # t c , which gives 

P [H t ] < P[K c t U H t ] < P[A7] + P[H t ]. 

By [201 Theorems 1 and 2] and by taking the union bound over all time steps, we have 

P[Kf] < texp(-ciL). (24) 



We will now derive an upper bound for P[flj]. Let £ > be a sufficiently small constant such that 
(1 — £)A > A c . Take the cube Q^l and tessellate it into cells of side length £, where I = C\ log 3//rf i, 
for C\ a sufficiently large constant in order to satisfy the assumptions of Proposition 15.21 Call a 
cell dense if it contains more than (1 — £)\£ d nodes. For 5 > 0, let D be the event that all cells 
inside Q2L are dense for at least (1 — 5)t time steps. Applying Proposition 15.21 we obtain a constant 
C2 such that 

P [D] > 1 - exp (-C2 



■ lQg 3+6/d t 

We use the event D to obtain an upper bound for P[-Ht] via 

P[H t ] < P[H t n D] + P[D C ] < P[H t DD]+ exp f-c 2 - ) . 



'log3+6A* t 

On the event D, by definition, we can find a collection S of (1 — <5)i time steps for which all cells of 
side length t are dense inside the cube Q2L- We set A = C 2 £ 2 for some sufficiently large constant 
C2. We define n as the first time step for which all cells of Q2L are dense. We now define Ti+i 
recursively as the first time step after Tj + A for which all cells are dense. Obviously, t\ < r 2 < • • • 
and if we take k = c^t/A = c' 3 t/ log 6 ^ d t for some constant C3, then we can ensure that on D we 
have Tfc < t — 1. 

For each i, let Ai be the event that u does no£ belong to a crossing component of at time Tj + A. 
Since when D holds we have Tfc < t — 1, we can write 



P[H t nD] < P 



f| Ai n D 



Li=l 



(26) 



For each i, let be the cr-field induced by the locations of the nodes of Ho from time to Tj. We 
now claim that for t sufficiently large there exists a positive constant C4 such that 

P [A4 I Fi] < e~ CA . (27) 
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We will define two events E\ , E2 such that for any F G T% we have 

P [Ai I F] < P [£? I F] + P [E§ I F] . 



(28) 



Take e > sufficiently small so that (1 — e)(l — £)A > A c , and let 5 be an independent Poisson 
point process of intensity (1 — e)(l — £)A. We define the events 

E\ = {u belongs to a crossing component of S in Ql} and 
E2 = {3 a coupling of H and n ri+ A so that S C n Ti+ A in Ql}, 

where "5 C n ri +A in Ql" means that the nodes of 5 that lie inside the cube Ql are a subset of 
the nodes of II Ti +A- 

Note that when E\ and E2 both hold, then u belongs to a crossing component of Ql at time Tj + A, 
which implies that A% does not hold. Since the intensity of S is strictly larger than A c and E\ is 
independent of F by construction, we obtain P \E\ \ F] > C5 for some constant C5 G (0, 1) by [20} 
Theorem 1]. 

All cells are dense at time Tj, by the definition of t, . Taking K and iT' appearing in Proposition 15. II 
to be K = 2L and K' = L, we see by the choice of A that for large enough t the condition for K' 
in Proposition 15.11 is satisfied and thus we obtain, uniformly over all F G JFi, that, for a positive 
constant Cq, 

P[^|F]<exp(-c 6 Alog 3 t). 
Plugging everything into (|28p we get 

P [A I F] < 1 - c 5 + exp (-c 6 A log 3 1) , 

which can be made strictly smaller than 1 by taking t sufficiently large. This establishes (1271) . 

Note that by definition we have Tj + A < t,i + i for all i, which gives Ai G J-i+i- We can write (I26p 
as 



-4; 



i=2 



u • 



P[i? ( nD] < P 

which by ([2?]) translates to 

P[f/t nD]< exp {-c A k) < exp ^ 
for a positive constant C7. Plugging this into (|25p concludes the proof of Theorem 11.61 



-c 7 - 



□ 



5.1 Coupling 

In this section we give the proof of Proposition 15.11 We begin by stating and proving a small 
technical lemma that will be used in the proof. 



Lemma 5.3. Assume e G (0, 1) and p > 0. Let A > 16d 2 p 2 /e 2 and R > 2^dA log^cfe" 1 ). Define 

5(z)= (27Aj^ eXP ( 2A— 

on Then we have 

/ 3(2) dz > 1 - e/2. 
Jb{o,r) 
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Proof. Let = exp (-^^) , for x G R. 

Note that Et=i(NI + /°) 2 = Hi + 2 P\\ z \\i + P 2 d > (\\z\\<i + pf , so 

d 

Y[i>(zi)<g(z), for z = ( Zl ,...,z d )£R d . (29) 
i=l 

Next observe that 

£^ = i-£^ <Bp (4)*>i- 7 g s >i-^>i-i i . 

By the Gaussian tail bound we have that 

^^)^<exp^j<^<-, 

for any e G (0,1)- Thus f*£f*-ij;(x) dx > 1 - ^j. Since [-R/Vd,R/Vd] d C B(0,R), we deduce 
from ([29]) that 

/ g(z)dz> [ T\^{zi)dz> (l-^) d > l-e/2. 

JB(0,B) J[-i?/v / d,i?/\/rfl' 1 ~ = i V 

□ 



We now proceed to the proof of Proposition 15.11 



Proof of Proposition \5.1[ We will construct H via three Poisson point processes. We start by 
defining Ho as a Poisson point process over Qk with intensity (1 — e/2)/3. Recall that $o has at least 
(3£ d nodes in each cell of Qk- Then, in any fixed cell, Ho has fewer nodes than «3>o if "0 has less than 
(3£ d nodes in that cell, which by a standard Chernoff bound (cf. Lemma [A. 1|) occurs with probability 
larger than 1 - exp ^_ £/2 (W 2 )^ d (! _ e '/ 3 ^ f or e ' such that (1 + e')(l - e/2) = 1. Since e G (0, 1) 

we have e' G (e/2, 1), and the probability above can be bounded below by 1 — exp (— ce 2 (3£ d ^ for 
some constant c = c{d). Let {Ho d> ^o} be the event that Ho has fewer nodes than <l?o in every cell 
of Qk- Using the union bound over cells we obtain 

P [H < *o] > 1 " exp(-ce 2 /3^). (30) 

If {Ho if* $o} holds, then we can map each node of Ho to a unique node of $o i n the same cell. We 
will now show that we can couple the motion of the nodes in Ho with the motion of their respective 
pairs in $o so that the probability that an arbitrary pair is at the same location at time A is 
sufficiently large. 

To describe the coupling, let v' be a node from Ho located at y' G Qk, and let v be the pair of v' 
in $o- Let y be the location of v in Qk, and note that since v and v' belong to the same cell we 
have \\y — y'\\2 < \fdl- We will construct a function g(z) that is smaller than the densities for the 
motions of v and v ! to the location y' + z, uniformly for z G K d . That is, 

1 ( maxill^lln, HV + z — •ullolA . s 

9(2,5 (2^P eXP ( "J < 31 » 
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= (2^7* GXP I 2A J • (32) 



for all z G Mr. 
To this end we set 



Note that this definition satisfies (|3ip since by the triangle inequality + z — 2/H2 < IW ~ ylb + ll^lb 
and \\y' — y\\2 < Define if) = 1 — J Rd g(z) dz. Then, with probability 1 — tp we can use the 

density function to sample a single location for the position of both v and v' at time A, and 
then set E' to be the Poisson point process with intensity (1 — ip)(l — e/2)/3 obtained by thinning Ho 
(i.e., deleting each node of Ho with probability ip). At this step we have crucially used the fact that 
the function g(z) in (|32p is oblivious of the location of v and, consequently, is independent of the 
point process $o- (If ° ne were to use the maximal coupling suggested by (I3ip . then the thinning 
probability would depend on $0, an d H would not be a Poisson point process.) 

Let H A be obtained from H after the nodes have moved according to the density function jr^p. 
Thus we are assured that the nodes of the Poisson point process H' A are a subset of the nodes of 
$A and are independent of the nodes of $01 where <£>a is obtained by letting the nodes of $0 move 
from time to time A. 

By Lemma 15.31 we get that if A and K — K' are large enough, then the integral of g{z) inside the 
ball B = -B(0, (K — K')/2) is larger than 1 — e/2. (We are interested in the ball B since for all 
z G Qk' we have z + B c Qk-) 

When {Ho if* $0} holds, E A consists of a subset of the nodes of Note that E A is a non- 
homogeneous Poisson point process over Qk- It remains to show that the intensity of H A is strictly 
larger than (1 — e)/3 in Qx> so that S can be obtained from H A via thinning; since H A is independent 
of <I>0) so is H. 

For z G M d , let n(z) be the intensity of H' A . Since H has no node outside Qk, we obtain for any 

z G 

M*) > (1 - V)(l - e/2)/3 / dx = (l- e/2)/3 / 5 (x) dx, 

where the inequality follows since z + B C for all z G From Lemma 15. 3| choosing 

the constants c\ and C2 sufficiently large we have j B g(x)dx > 1 — e/2. We then obtain ^(2) > 
(1 — e/2) 2 /3 > (1 — e)/3, which is the intensity of E. Therefore, when {Ho ^ $0} holds, which occurs 
with probability given by (I30p . the nodes of H are a subset of the nodes of $a, which completes 
the proof of Proposition 15.11 □ 



5.2 Density 

In this section we prove Proposition 15.21 using a multi-scale argument. Since the argument is rather 
involved, we begin with a high-level overview. 

Proof overview 

Our goal is to show that if we tessellate the cube Ql, with L = t, into cells of volume of order 
(logt) c , then the probability that all cells contain sufficiently many nodes for a fraction 1 — e of the 
time steps is at least the expression given in Proposition 15.21 

We start at scale 1 with the cube Ql x where L\ > L. We tessellate Ql 1 into cells that are so large 
that we can easily show that with very high probability during all time steps all these cells contain 
sufficiently many nodes. We refer to this as the event that "the density condition is satisfied at 
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1 h-, <M 



(a) 



(b) 



(c) 



Figure 1: (a) The cube Ql _ 1 and the smaller cube Q/, . obtained when going from scale j — 1 to 
scale j. (b) The tessellation of into cells of side length £7-1. (c) The finer tessellation of Ql^ 

into cells of side length tj < 

all steps for scale 1." Then, when going from scale j — 1 to scale j, we take a smaller cube Qx,. 
with < Lj—i, and tessellate it into cells that are smaller than the cells at the previous scale (see 
Figure [TJ . We define the density condition for scale j at a given time step as the event that all the 
cells at scale j contain a number of nodes that is sufficiently large but strictly smaller than the one 
used for the density condition for scale j — 1. Since this density requirement becomes less strict 
when going from scale j — 1 to scale j, we will be able to show that the density condition for scale j 
is satisfied for a large fraction of the time steps at which the density condition is satisfied for scale 
j — 1. We repeat this procedure until we obtain, at the last scale, the cube Ql and cells of side 
length I. 

The importance of the multi-scale approach is that it allows us to recover quickly from instances 
of low density, i.e., if the density condition holds in scale j — 1 but fails (at some time) in scale j, 
there are enough nodes nearby to recover density shortly thereafter. 

We now proceed to the detailed argument. 
Full proof 

Let k be the number of scales; we will see in a moment that k = O(logi) will suffice. Let L\ > 
L2 > • • ■ > L K = L such that L\ = t 2 and L K = t. 

Let i\ > £2 > • • • > = £■ At scale j, we consider the cube Ql^ and tessellate it into cells of 
side length £j (see Figure QJb-c) ) . We say that a cell is dense at a given time step for scale j if it 
contains more than (1 — nodes at that step, where the £j satisfy 

2 = 6 < 6 < ■ ■ ■ < 6t = £ and - = _ for all j. 

We start by analyzing the event that all cells are dense for scale 1 during all time steps, which we 
denote by D\. The next lemma shows that D\ occurs with very high probability. 

Lemma 5.4. If if > Clogt for some large enough constant C, then there exists a constant c such 
that 

P [£>i] > 1 - exp f-cXifj . 

Proof. For any fixed time i and cell k, the number of nodes in k at time i is given by a Poisson 
random variable with mean Xif. Then, using a Chernoff bound (cf. Lemma lA.lj) . we obtain that 
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Figure 2: Illustration for how 1 time interval of scale j — 1 gives 4 subintervals of scale j. 

there are more than (1 — £i)\£f nodes in that cell at that time step with probability larger than 
1 - exp(-^A£f/2). The number of cells inside Q Ll is 0(t 2d ) by our choice of L\ and £\. The proof 
is completed by taking the union bound over all cells and time steps, and using the assumption on 
£ v □ 

We will need to disregard some time steps when going from one scale to the next. During this 
discussion it will be useful to refer to Figure [2j Let Sj be the number of time steps considered 
for scale j. We start with si = t so that at scale 1 all time steps are considered; we will have 
s\ > S2 > • • • > s K . For each scale j, we will split time into intervals of rrij consecutive time steps. 
We start with mi = t, so that at scale 1 we have only one time interval of length t. 

In each interval [a, a + rrij-i) at scale j — 1 we consider the following four separated subintervals of 
length rrij (see Figure [2]): 

[a + kAj_i + (fc — 1)7717, a + kAj_i + krrij), for k = 1,2,3,4, (33) 

where 

rrij = J —. (34) 

We will set the Aj in a moment. We skip A,_i steps in order to allow the nodes to move far 
enough and enable the application of the coupling from Proposition 15.11 Note that this gives 

For a given scale j, we say that a time interval is dense if all cells are dense during all the time 
steps contained in this time interval, i.e., each cell contains more than (1 — £j)\£j nodes at all time 
steps. 

Let = ei < £2 < • • • < e K = e satisfy €j — ej-\ = For each scale j > 1, we define the event 

Dj = {a fraction of at least > ^1 — time intervals of scale j are dense}. (35) 

If D K holds, the number of time steps for which all cells are dense for the last scale k is at least 

-f)-(-l)-n(-^)^(-l)'(-|^)- <»> 



Since we are aiming to obtain (1 — e)t time steps for which the density condition is satisfied for the 
last scale, we set Aj to satisfy 

- = r < 37 > 
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for all j. The value of Aj must be sufficiently large to allow nodes to move over a distance £j. We 
then define £j by 

Aj = C'£ 2 k 2 , (38) 

where C is a sufhciently large constant. 
From ((Ml), (HU) and ((38]), we obtain 



8C'k 3 3+1 3 V 4 8k 



and £ j+1 = £ n ll-± (39) 



Since mi = t, we get that £\ = 8C ^ k3 f < g^ri and since we want to get £ d K = £ d > C(logi) 3 , it is 
easy to see that k = O(logt) is sufficient. 

For any time step i, let T% be the cr-field induced by the locations of the nodes of Ho from time 
up to time i. 

Lemma 5.5. Let A = [a, a + rrij) be a time interval considered in scale j. We write b = a — Aj_i 
and E = {at time b all cells are dense for scale j — 1}. Let £ d > C(logt) 3 for some sufficiently 
large constant C > 0. Then there exists a constant c such that 

P [A not dense, E | J-j,] < exp (-c\£ d /K 2 

Proof. For any F £ Jj such that F n E = the lemma clearly holds. We then take F n E ^ 
and give an upper bound for P [A not dense | E, F] . Let $6 be the point process obtained at time 
6 after conditioning on F n E. We first fix a time w G A and derive an upper bound for 

P [at time w not all cells are dense for scale j \ E,F] . 

Since we condition on E, all cells are dense for scale j — 1 at time b. We now set 5 such that 
(1 — S) 2 (l — = 1 — £j, which implies 5 = — Cj-i)- We also choose a constant c and the 

constant C appearing in the definition of A,- in (|38|) so that, setting 



Lj < Lj-i - cy Aj_i log -, (40) 

allows us to apply Proposition 15.11 with K = Lj^i and K' = Lj. Thus we obtain a fresh Poisson 
point process 3 with intensity (1 — S)(l — £j-i)A that can be coupled with & w (which is the point 
process obtained at time w after the points of have moved for time w — b) in such a way that S 
is stochastically dominated by Q w inside Qlj with probability at least 

1 _ exp (_ Cl 5 2 (l - C-OA^i) , (41) 

for some positive constant c\. We note that the choice of L\ = t 2 and the fact that k = O(logt) 
together with equation f)38|) gives that it is always possible to choose the Lj's satisfying (I40p and 
such that L K = t. 

A given cell is dense for scale j at time w if S contains at least (1 — £j)\£j nodes in that cell, which 
by the choice of S happens with probability at least 1 — exp ^— C2^ 2 (l — S)(l — £j_i)A£^ for some 

constant C2 (cf. Lemma |A. 1|) . The proof is completed by taking the union bound over all cells and 
over all time steps in A and using the condition for £. □ 
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We now use Lemma 15.51 to give an upper bound for P 
was denned in (1351). 



that holds for all j, where Dj 



Lemma 5.6. If l d > C(logi) 3 for some large enough C , then there exists a constant c such that 
for any j > 2 we have 



Proof. If Dj-i happens, then there are at least (l — -^p) ^. 1 dense time intervals for scale j — 1. 
When we go to scale j, these intervals will give us 



4 1 



2 / m 



(42) 



time intervals that we will consider for scale j. On the other hand, if the event D c - holds, then 
there are less than 



2 / mi 



(43) 



dense intervals for scale j. Let w be obtained by subtracting 



from 



that is, 



w 



£ 7 - e 7-l 



Let Z be the number of subintervals [a, a + nij) of scale j that are not dense for scale j, but are 
such that the time step a — A,_i is dense for scale j — 1. (We call a time step dense if all cells are 
dense at that time.) It is easy to see that if both and D c - happen, then Z > w. 

We can write Z as a sum of Sj/rrij indicator random variables one for each time interval of scale 
j. Although the l^'s depend on one another, Lemma [5.51 gives that the probability that Ik = 1 given 
an arbitrary realization of the previous k — 1 indicators is smaller than pj = exp ^— ci£ 2 A^/k 2 ^ 

for some constant ci. Therefore, Z is stochastically dominated by a random variable Z' obtained 
as a sum of sj/rrij i.i.d. Bernoulli random variables with mean pj. Using a Chernoff bound (cf. 
Lemma lA.2j) . we obtain 



P [Z 1 



> w 



Z'-B \Z'] > ^- 



< 



exp 



s i / e i - i 



mi 



log 



e 7 _i 



2pi 



(44) 



Note that ej - ej-i = ^ and log^J 1 ) = G(^ 2 A^/k 2 ). Also £<j > l d > C(logt) 3 and k = O(logt) 
so we obtain a constant C2 such that 



P \Z' >w]< exp -c 2 ^ 2 As : 



mi k° 



Recall from ([39]) that ^- = s q, k 3 ■ By ([36|) and ([37|) we have that = Q(t) for all j, so we finally 
obtain 



P [Z' >w]< exp 



2f2 



-C 3 - 



■\£ d - 2 t 



for some constant C3 > 0. Using £j > I and the assumption on I in the statement of the lemma 
completes the proof. □ 
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We are now in a position to prove Proposition 15.21 



Proof of Proposition \5.2[ To prove Proposition 15.21 we need to derive an upper bound for 
P [D%\. Note that P [D%] < P [D% n Ae-i] + P • Applying this inequality recursively for the 

term P [-D^J we obtain 

K 

Each term in the sum can be bounded using Lemma 15.61 and the last term can be bounded using 
Lemma 15.41 The proof is completed since k = O(logt) and the initial value 

□ 



6 Broadcast time 

In this section we use Theorem 11.61 to prove Corollary 11.71 for a finite mobile network of volume 
n/A. 

We may relate the mobile geometric graph model on the torus to a model on R rf as follows. Let S n 
denote the cube Qt n /x\V<i- The initial distribution of the nodes is a Poisson point process over R rf 
with intensity A on S n and zero elsewhere. We allow the nodes to move according to Brownian 
motion over M. as usual, and at each time step we project the location of each node onto S n so 
that nodes "wrap around" S n when they reach the boundary. 

Proof of Corollary \1.7[ Let t = Clogn(loglogra) 3+6 / d for some sufficiently large constant C = 
C{d). We define a giant component as a connected component that contains at least two nodes 
at distance larger than ( W//- ^ — . It follows from |2CH Theorem 2] and the union bound over time 
steps that, with probability 1 — e -®( nl/d ) ) Q i contains a unique giant component for all integer 
i e [0,2i-l]. 

The proof proceeds in two stages. First, we show that for any fixed i £ [0, 24 — 1], w.h.p. the giant 
component of Gi has at least one node in common with the giant component of Gj+i. This means 
that, once the message has reached the giant component, it will reach any node v as soon as v itself 
belongs to the giant component. Then we show that, after 4 steps, all nodes have belonged to the 
giant component w.h.p. This implies that broadcast is achieved after 24 steps w.h.p. 

To establish the first stage, let e > be sufficiently small so that (1 — e)A > A c . We use the thinning 
property to split H into two Poisson point processes, IL^ and IT/, with intensities (1 — e)A and eA 
respectively. Let G\ and G' i+1 be the graphs induced by Ll'j and LL[ +1 respectively. Then with 
probability 1 — e -®( nl/d ) both G\ and G' i+1 contain a unique giant component \20\ Theorem 2]. We 
show that at least one node from LI" belongs to both giant components. For any node v of LI'/, 
the probability that v belongs to the giant component of G\ is larger than some constant c = c(d). 
Moreover, using the FKG inequality we can show that v belongs to the giant components of both 
G\ and G' i+l with probability larger than c 2 . Therefore, using the thinning property again, we can 
show that the nodes from LI" that belong to the giant components of both G' i and G' i+1 form a 
Poisson point process with intensity eAc 2 , since c does not depend on IT/. Hence, there will be at 
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least one such node inside S n with probability 1 — e ec n , and this stage is concluded by taking the 
union bound over time steps i. 

We now proceed to the second stage of the proof. We first need to show that the tail bound on 
T perc from Theorem 11.61 also holds when applied to the finite region S n defined above. Note that 
all the derivations in the proof of Theorem 11.61 were restricted to the cube Ql^ where L\ = t 2 was 
defined in Section 15.21 We have that Q l ± is contained inside S n for all sufficiently large n since 
L\ = t 2 = (9(log 2 n(log log n) 6+12 / rf ) while S n has side length (n/X) l l d . In order to check that the 
toroidal boundary conditions do not affect the result, it suffices to observe that, during the time 
interval [0, 2t], no node moved distance larger than ( n /^ — w.h.p. 

Now note that, by a Chernoff bound, G has at most (1 + S)n nodes with probability larger than 
1 — e~ n ( n ^ for any fixed 5 > 0. These nodes are indistinguishable, so letting p be the probability that 
an arbitrary node has percolation time at least t, we can use the union bound to deduce that this 
applies to at least one node in G with probability at most (l + d)np. Let v be an arbitrary node. In 
order to relate p to the result of Theorem II. 6^ we can use translation invariance and assume that v 
is at the origin. Then, by the "Palm theory" of Poisson point processes [23], p is equivalent to the 
tail of the percolation time for a node added at the origin, which is precisely P [T pcrc > t]. Thus 
finally, using Theorem 11.61 we get p < exp(— c ( logt *3+6/d ); which can be made o(l/n) by setting C 
sufficiently large in the definition of t. 

We then obtain that with probability 1 — o(l/n) all nodes of G have been in the giant component 
during the time interval [0,t — 1], which implies that at time step t — 1, the nodes of the giant 
component contain the message being broadcast. By stationarity, with probability 1 — o(l/n) all 
nodes have been in the giant component during the time interval [t, 2t — 1], and thus have received 
the message by time It. This completes the proof of Corollary 11.71 □ 

Remark 6.1. It is easy to see that the above result also holds in the case where the graph 
has exactly n nodes. The proof above shows that, by setting C large enough, we can ensure 
P P~bc > 2t] = o(l/n) for the given value of t. Also, it is well known that a Poisson random 
variable with mean n takes the value n with probability p = Q(l/ \/n). Therefore, for a graph with 
exactly n nodes, we have Pr[Tt, c < t] = p °^ n ^ = 1 — o{l/^/n). 
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A Standard large deviation results 

We use the following standard Chernoff bounds and large deviation results. 

Lemma A.l (Chernoff bound for Poisson). Let P be a Poisson random variable with mean A. 
Then, for any < e < 1, 



Lemma A. 2 (Chernoff bound for binomial [1] Corollary A. 1.10]). Let X be the sum of n i.i.d. 



Bernoulli random variables with mean p. Then, P [X > np + a] < exp I a — (pn + a) log ( 1 + — 




and 



P [P < (1 - e)A] < exp 
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